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Abstract 
rl , 

'■^ The presently observed cosmological baryon asymmetry has been finahy de- 

termined at the time of the electroweak phase transition, when baryon and 
^ ■ lepton number violating interactions fell out of thermal equilibrium. We dis- 

cuss the thermodynamics of the phase transition based on the free energy of 
the SU(2) Higgs model at finite temperature, which has been studied in per- 
turbation theory and lattice simulations. The results suggest that the baryon 
asymmetry has been generated by lepton number violating interactions in the 
symmetric phase of the standard model, i.e., at temperatures above the critical 
temperature of the electroweak transition. The observed value of the baryon 
asymmetry, n^/s ~ 10"^'', is naturally obtained in an extension of the standard 
model with right-handed neutrinos where B — L \s broken at the unification 
scale Aqut ~ 10^^ GeV. The corresponding pattern of masses and mixings of 
the light neutrinos z^e, and v-r is briefly described. 



*To be published in Quarks '96, Proc. of the IXth International Seminar, Yaroslavl, Russia, 1996 



1 Introduction 



In the standard model of electroweak interactions all masses are generated by the 
Higgs mechanism. As a consequence, at high temperatures a transition occurs from 
a massive low-temperature phase to a 'massless' high-temperature phase, where the 
Higgs vacuum expectation value 'evaporates' and the electroweak symmetry is 're- 
stored' ig. 

Due to the chiral nature of the weak interactions baryon number (B) and lepton 
number (L) are not conserved in the standard model 0. At zero temperature this 
has no observable effect due to the smallness of the weak coupling. However, as 
the temperature approaches the critical temperature of the electroweak phase 
transition, B and L violating processes come into thermal equilibrium |^ . Their rate 
is determined by the free energy of sphaleron-type field configurations which carry 
topological charge. In the standard model they induce an effective interaction of all 
left-handed fermions (cf. Fig. 1) which violates baryon and lepton number by three 
units, 

AB = AL = 3. (1) 
Since B and L violating processes fall out of thermal equilibrium below T^, the 




Figure 1: One of the 12-fermion processes which are in thermal equilibrium in the 
high-temperature phase of the standard model. 



presently observed value of the baryon asymmetry of the universe has finally been 
determined at the electroweak transition. Hence, the study of the thermodynamics 
of this transition is of great cosmological significance. 

Sphaleron processes conserve B — L. In the high-temperature, 'symmetric' phase, 
where B and L violating processes are expected to stay in thermal equilibrium over 
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a wide range of temperatures, the expectation values of B and L are related, 



{B)t^C{B-L)t^^^{L)t, (2) 

where C = || in the standard model 0. Hence, baryon and lepton number are 
correlated in the symmetric phase, and the generation of a baryon asymmetry requires 
lepton number violation. 




Figure 2: Free energy as function of the order parameter at a temperature slightly 
above the critical temperature of a first-order phase transition. 



In the abelian Higgs model, i.e., scalar electrodynamics, the high-temperature 
phase transition is known to be of first order for sufficiently small Higgs masses ^. 
The expectation value of the Higgs field plays the role of an order parameter. At 
temperatures just above the critical temperature the Higgs phase is metastable. It 
decays in a first-order transition to the symmetric phase, which is accompanied by a 
jump in the order parameter Av = Vb — Vg (cf. Fig. 2). 

A similar behaviour is expected in the standard model, with a critical temperature 

— 1/2 

given by the Fermi constant, Tc ~ Gp . In a first-order phase transition a strong 
deviation from thermal equilibrium can occur. Since the standard model contains 
CP violating and, at high temperatures, also B violating interactions, all conditions 
for baryogenesis are fulfilled. It is therefore conceivable that the cosmological baryon 
asymmetry has indeed been generated at the electroweak phase transition 0. Such 
a scenario requires that a produced baryon asymmetry is not erased by the sphaleron 
processes in the Higgs phase close to the critical temperature. From this condition a 
lower bound on the jump of order parameter at the phase transition can be derived 
(see 0), 
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> 1.2. 



(3) 



This condition is necessary for electroweak baryogenesis although far from sufficient. 
A discussion of the comphcated nonequihbrium processes in the electroweak plasma 
can be found in [0. 

In order to examine whether condition (j^) is satisfied, one has to study the ther- 
modynamics of the electroweak transition near the critical temperature. Here the 
main obstacle are the well known infrared divergencies of non-abelian gauge theo- 
ries at high temperature. As a first step towards a treatment of the full standard 
model the following simplifications are usually made: fermions are integrated out 
using high-temperature perturbation theory, and the electromagnetic interaction is 
neglected (sinB^y = 0). One is then left with the SU(2) Higgs model. 



2 Thermodynamics of the SU(2) Higgs model 



The observables which characterize a first-order phase transition are critical tem- 
perature, jump in the order parameter and latent heat. Also important are surface 
tension and correlation lengths which are more difficult to compute and which we 
shall not discuss in detail. Our analysis will be based on the gauge invariant 'order 
parameter' ($^^$) and the corresponding free energy computed in lattice simulations 
and perturbation theory. In our discussion we shall closely follow [Q. 

It is well known that the electroweak phase transition is infiuenced by non- 
perturbative effects whose size is governed by the confinement scale of the effective 
three-dimensional theory which describes the high-temperature limit of the SU(2) 
Higgs model. These effects are particularly relevant in the symmetric phase, and one 
may worry to what extent a perturbative analysis of the phase transition can yield 
sensible results. As we shall see, perturbation theory is self-consistent at two-loop 
order. A comparison with results obtained by lattice simulations will then enable us 
to estimate the effect of non-perturbative corrections. 

Generalities 

The action of the SU(2) Higgs model at finite temperature T reads 



Sp[^,W]= [ 

J 3 




(4) 



with 




1 



1 



(5) 



(6) 
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Here PF^ is the vector field, a is the Higgs field, vr is the Goldstone field, f is the 
triplet of Pauli matrices, and VL is the spatial volume. For perturbative calculations 
gauge fixing and ghost terms have to be added to the action (^). 

The free energy density of the system, W{T, J), is given by the partition function, 
i.e., the trace of the density matrix. 



exp{-pnW{T, J)) = Tr exp 



-(3 {h + J j^d^x^^^ 



(7) 



where H is the Hamilton operator of the theory, and $ is the operator describing the 
Higgs field. We have added a source J, with 5^ J = 0, coupled to the spatial average 
of the gauge invariant composite operator <l>^<l> (here and below the trace operator 
acting on $^<l> is omitted for brevity). The partition function can be expressed as a 
euclidian functional integral (see 0]), 

exp (-/5fiVr(T, J)) = j D^D<^^DW^exp(^- J ds(L + J<l>"f$)^ , (8) 

where L is the euclidean lagrangian density, and the bosonic fields $ and Wf^ satisfy 
periodic boundary conditions at r = and t = [3. Eq. (||) is the starting point of 
perturbative as well as numerical evaluations of the free energy. 

Note, that the source J in Eq. (|]) couples to a gauge invariant composite field. 
Hence, the free energy W{T, J) is gauge independent. The spatially constant source 
J simply redefines the mass term in the action (^. This is in contrast to the usually 
considered generating function of connected Green functions at zero momentum, 



exp (-pVtWiT,]- J)) = J^D^D^^DWf^exp j^dx (l + J$^$ + ja) 



(9) 



Here the source j couples to a gauge dependent quantity, the field a. Consequently, 
W{T,i]J) is gauge dependent and not a physical observable. We have also kept 
the dependence on the source J. From W(T,j; J) one obtains the effective potential 
V{T,(f; J) via Legendre transformation, with (p = dW/dj. The wanted free energy 
density W{T, J) can now be obtained from the effective potential V. In the infinite 
volume limit, one has 

W{T, J) = V{T, (^™„(T, J), J) , (10) 

where (frnin{T,J) is the global minimum of the effective potential V{T,(f;J). For 
arbitrary values of ip the potential V is gauge dependent. However, its value at the 
minimum is known to be gauge independent, yielding a gauge independent free energy 
W{T, J). 

At the critical temperature of a first-order transition the order parameter p, 
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and the energy density E, 



E{T, J) = W{T, J) - T-^W{T, J) , (12) 

are discontinuous. The jump in the energy density is the latent heat AQ. 

For the first-order phase transition from liquid to vapour there exists a well known 
relation between the latent heat and the change of the molar volume, the Clausius- 
Clapeyron equation In the electroweak phase transition the order parameter 



($^$) plays the role of the molar volume, and a completely analogous relation can 
be derived. 

The electroweak plasma can exist in two phases, the massive low-temperature 
Higgs phase with free energy Wb{T, J) and the massless high-temperature symmetric 
phase with free energy Ws{T, J). In the J — T-plane the boundary between the two 
phases is determined by the equilibrium condition 

Ws{T,J{T)) = Wh{T,J{T)) . (13) 

This equilibrium condition yields a useful connection between the latent heat AQ and 
the jump in the order parameter Ap. Using the definitions 

AQ = -T— {Ws -W,) , Ap = 2 ^ {Ws - W,) , (14) 
one easily obtains 

AO = lApr;|^. (15) 

This is the Clausius-Clapeyron equation of the electroweak phase transition. Together 
with a dimensional analysis it implies 

AQ = -]^mlAp{l + 0{g\\)) , (16) 

where niH = v^— 2/?(l + 0{g'^,\)) is the physical Higgs mass at zero temperature. 
This relation provides a useful check for perturbative as well as lattice results. 

Perturbation theory 

Near the ground state, J = 0, the free energy W{T, J) can be evaluated as power 
series in the couplings g and A by means of resummed perturbation theory which has 
been carried out up to two loops [|ll|, |12|, |[. Here, thermal corrections are added to 



the tree-level masses of the scalar fields and the longitudinal component of the vector 
boson field, 

5Sp = pj d'x (iaoiT2(a2 + n') + ^a^T'Wl^ . (17) 
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The sum of tree-level masses and thermal corrections then enters the boson propa- 
gators in loop diagrams, and 5Sp = —6S13 is treated as counter term. It turns out 



that the resummation of static modes only is a preferred procedure [11|. Hence, in 
Eq. (|T7|) the fields a and vr do not depend on the imaginary time r. To leading order 
in the couplings, one obtains for the parameters in Eq. ( ]T7| ) from one-loop self energy 
corrections Oqi = j^g'^ + |A, ai = ^g"^. 

The masses of the boson propagators are obtained from Eqs. (||) and by 
shifting the Higgs field a by the average field ip. This yields mi, rriT, and m.,^ for 
longitudinal and transverse part of the vector field, the Higgs field and the Goldstone 
field, respectively. The resummation procedure can be optimized by adding terms of 
higher order in the couplings to Eq. (p^. In the Higgs phase, where one is only inter- 
ested in the effective potential close to the minimum, the choice of a field dependent 
correction 55*^ is useful which yields for the scalar masses. 



mi = 2Xifi' , mi = . (18) 

Hence, no thermal resummation for scalar masses is performed in the Higgs phase. 
In the symmetric phase it is useful to determine the scalar masses self-consistently by 

ml = ml = ^-^V{T,0;J) . (19) 

For given vector boson masses m^ and m-p, this is a gap equation for the scalar 
masses, which can be solved at each order of the loop expansion. 

Having specified the vector boson and scalar masses in the Higgs phase and in 
the symmetric phase, the free energy can be calculated from the effective potential 
V, using the relation (p!0|). Its two local minima yield the free energy Ws{T, J) and 
Wi,{T, J) in the symmetric phase and the Higgs phase, respectively. The free energy 
of the ground state is 

W{T, 0) = min{iy,(T, 0), W,{T, 0)} . (20) 

It is a concave function whose derivative is discontinuous at the critical temperature 
Tc. We can also consider the dependence of the free energy on the external source 
at the critical temperature Tc. This function, W{Tc, J), is shown in Fig. 3. Here we 
have subtracted in both phases the huge linear term T'^J/6, whose sole effect is to 
shift the expectation value ($'''$) by 

From the free energy W{T, J) one can obtain the gauge invariant effective potential 
V{T,p) by means of a Legendre transformation. Since the derivative of W{T, J) is 
not continuous everywhere, one has to use the definition (see 



V{T,p) = snp{W{T,J)-lpJ} . (21) 
J 2 
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Figure 3: Free energy at the critical temperature as a function of the source J {niH = 70 GeV). 
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Figure 4: Gauge invariant effective potential (solid line) and its analytic continuations from the 
single phase regions into the mixed phase region (dashed line). The cross denotes the matching 
point. 
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Figure 5: Jump of the order parameter ipi- 



/ Ap in units of the critical temperature T^. 



This yields the convex, non-analytic function, which is plotted in Fig. 4 as full line. 
One may also compute the ordinary Legendre transform Vs(T,p) and Vb(T,p) of 
Ws(T, J) and Wb(T, J), respectively. Vs and Vb are also shown in Fig. 4. In the region 
outside of the two local minima, Vs and Vb, respectively, coincide with the convex 
effective potential V{T,p). Between the two local minima, Vg and Vb represent two 
analytical continuations of V(T,p), which meet at the 'matching point' pm = 
At this point, marked by a cross in the plot, the first derivatives of both curves 
coincide. 

The non-convex 'effective potential' obtained by combining Vg and Vb on both 
sides of the 'matching point' has a barrier between symmetric and Higgs phase like 
the ordinary effective potential V{T, ip; J). Note, however, the difference in the range 
of fields. For V one has < (f < oo, whereas for the potential V the field p varies in 
the range — oo < p < oo. The generation of a barrier between two local minima as 
analytic continuation from a convex effective potential is reminiscent of the treatment 
of first-order phase transitions in condensed matter physics . However, the precise 
physical meaning of the resulting non-convex 'effective potential' still remains to be 
understood. 



9 



Comparison with lattice simulations 



The asymmetric resummation described above has been carried out up to two- 
loop order. In the Higgs phase the scalar masses are given by Eq. (|T^) , and in the 
symmetric phase they are self-consistently determined from Eq. (p^Q]). Note, that in 
the two-loop calculation the counter term to be inserted in the one-loop graph is 
0{g^), whereas the scalar mass determined from Eq. (^) is of higher order in g. 

In the symmetric phase the self-consistently determined scalar masses are infrared 
divergent. With rriT = the two- loop potential yields a contribution which 

diverges logarithmically at ^ 0, 

~ - -^^T^ In fSruT . (22) 
1287r2 ^ ^ ' 



Following [|T^ one may regularize this divergence by means of a 'magnetic mass' term. 
In Eq. ( ^21) one substitutes m'^ = g'^(p'^/4 + '-f^g^T^/{9TT'^). In the following numerical 



results will be given for 7 = 1, which is obtained by one-loop gap equations |T^, |T( 
The results change only insignificantly if the parameter 7 is varied between 0.3 and 3.0. 
In addition to the resummation procedure one has to choose a renormalization scheme. 
A good choice is the MS-scheme with fi = T, supplemented by finite counter terms 
6^"^ and 6X which account for the most important zero-temperature renormalization 
effects. 

Given the two-loop potential in the symmetric phase and in the Higgs phase, one 
can numerically determine the critical temperature T^, where the two potentials at 
their respective local minima are degenerate. Differentiation with respect to tem- 
perature and the external source J then yields latent heat and jump in the order 
parameter p = 2($"''$). The result for (fc = V^P is shown in Fig. 5, labelled 'new 
resummation'. Here the zero-temperature standard model values = 80 GeV and 
(yf2 = 0.57 have been used. In the case of the 'old resummation' |]12| (f^ corresponds to 
the position of the second minimum. The one-loop potential is 0{g^, A^/^), the two- 
loop potential is 0{g'^,X'^). As Fig. 5 illustrates, the 'new resummation' procedure 
improves the convergence significantly. The relative change of an observable from 
one-loop to two-loop may be characterized hj 6 = 2\0i — 02\/{0i + O2). For the 
jump in the order parameter ipc, ^ increases from ~ 0.01 at rriH = 40 GeV to ~ 0.2 
at niH = 70 GeV. Above rriH ~ 80 GeV the convergence deteriorates rapidly, and the 
perturbative calculation is no longer self-consistent. 

Observables of the four- dimensional SU(2) Higgs model at finite temperature can 
be directly computed by means of Monte Carlo simulations on lattices with spatial 
size Lg and temporal size Lf -C Lg. Such finite temperature simulations at small 
values of A were initiated in [jl^, and high statistics simulations were performed in 
[IE, |TU[. For the two Higgs masses rriH ^ 18 GeV and ttih ~ 49 GeV the quantities 
Tc, AQ and ipc = vt were computed on = 2 and Lt = 3 lattices, and for tjih ~ 35 
GeV a detailed study of the first-order phase transition was carried out on lattices of 
temporal size up to Lj = 5. 
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In a complementary approach detailed studies have also been carried out based on 
dimensional reduction where non-zero Matsubara frequencies are first integrated out 
perturbatively. Recent results on dimensional reduction are described in pO, pT| where 



also references to previous work can be found. Numerical simulations were performed 
for the effective three dimensional theory which led to a quantitative description of 
the first-order phase transition for Higgs masses up to ~ 70 GeV p2|, E3 . 
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Figure 6: Jump of the order parameter at the critical temperature. Comparison of 
four-dimensional simulations (triangles, squares |T8|) with three-dimensional simula- 
tions (stars ||22|) and perturbation theory P]. From ||2^ 
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In Fig. 6 results for the jump in the order parameter, Av{Tc) = vt, are compared 
which were obtained by simulations of the four- dimensional theory , the three- 
dimensional theory |^2| and perturbation theory f^, respectively. The comparison is 
made for the parameter values Mw = 80 GeV and g = 0.57. The statistical errors 
of the numerical simulations are so small that they are invisible in the figure. The 
agreement between the three approaches is remarkable and certainly better than the 
systematic uncertainties of perturbation theory (cf. Fig. 5). The continuum limit for 
the critical temperature has been studied in detail in [|I^ for mn — 34 GeV , i.e., 
Rhw = Mu/My/ ~ 0.42. The results are shown in Fig. 7, where the error bars 
include the statistical error and an estimate of the systematic error. A comparison 
is made with perturbation theory and the three-dimensional simulation. The lattice 
simulations yield a critical temperature slightly below the result from perturbation 
theory, but the difference is not significant. 

The agreement between results from perturbation theory and non-perturbative 
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lattice simulations is surprizing, since in the symmetric phase perturbation theory 
is known to be infrared divergent, which prevents a straightforward extension of the 
present two-loop calculation to three loops. The infrared behaviour of a running 
gauge coupling has been studied in p6|. In the symmetric phase the scalar masses 
depend logarithmically on an infrared cutoff of order the magnetic scale ~ g^T. This 
is reminiscent of the Debye screening length in pure gauge theory which, at two-loop 
order, also requires an infrared cutoff ~ g^T . Here, a non-perturbative definition 
of the Debye screening length can be given such that the dominant contribution is 
given by the perturbative result and an additional non-perturbative contribution can 
be evaluated in a well-defined manner pSl. Such a split into a dominant perturbative 



contribution and a non-perturbative remainder should also be possible for the free 
energy in the symmetric phase. This would then justify the perturbative treatment of 
the first-order electroweak phase transition. For the pure gauge theory this problem 



has been discussed in 
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Figure 7: Numerical results for the ratio of critical temperature and Higgs boson mass 
versus {aTcY = L^^; the full line is the extrapolation to the continuum limit |T9|. The 



dashed horizontal line shows the prediction of two-loop perturbation theory p[. The 
hexagon between continuum limit and perturbation theory represents the result of a 
three-dimensional simulation P2[, adjusted to the parameters of the four- dimensional 



calculation 25 . 
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Perturbation theory in the symmetric phase fails with respect to correlation lengths 
of gauge-invariant vector and scalar fields. Lattice simulations yield shorter cor- 
relation lengths for vector fields in the symmetric phase than in the Higgs phase 
The contrary is true for the correlation length of the vector field 
a puzzle which still remains to be resolved. The 
latter screening lengths appear to be related to the critical Higgs mass ttih ~ 80 GeV 
where the first-order phase transition turns into a smooth crossover 



18, 22, m,m 23 



evaluated in a fixed gauge BTI, |32 



Independent of the intriguing non-perturbative features of the symmetric phase, 
which require further studies, it is now known that the necessary condition (|^) for 
electroweak baryogenesis is not satisfied in the standard model. For Higgs boson 
masses above the present experimental bound of 58 GeV one has according to 
Fig. 6, 



< 0.7 



(23) 



which is much smaller than the lower bound (^. For appropriate choices of parameters 
it is possible to satisfy this bound in some extensions of the standard model, but in 
these models it has not been demonstrated that electroweak baryogenesis is indeed 
possible. This strongly suggests that the baryon asymmetry has been generated in 
the high-temperature phase of the standard model, as it is the case for instance in 
conventional grand unified theories. 



3 Baryogenesis via leptogenesis 

In the high temperature phase of the standard model the asymmetries of baryon 
number and lepton number are proportional in thermal equilibrium (cf. (0)), 

{B)t ^ C{B - L)t ^ -^-j{L)t . 

In the standard model, as well as its unified extension based on the group SU(5), 
i? — L is conserved. Hence, no asymmetry in B — L can be generated, and (-B)t 
vanishes. Furthermore, as discussed above, baryogenesis at the electroweak phase 
transition appears unlikely. As a consequence, the non- vanishing of the baryon asym- 
metry is a strong argument for lepton number violation. This is naturally realized by 
adding right-handed Majorana neutrinos to the standard model. This extension of 
the standard model can be embedded into grand unified theories with gauge groups 
containing SO(IO) ||3^. Heavy right-handed Majorana neutrinos can also explain the 
smallness of the light neutrino masses via the see-saw mechanism |^ . 



The connection between baryon number and lepton number at high tempera- 
tures can be used to generate a baryon asymmetry. This was suggested by Fukugita 
and Yanagida [^. The primordial lepton asymmetry is generated by the out-of- 
equilibrium decay of heavy Majorana neutrinos in the standard manner. This mech- 
anism has subsequently been studied by several authors [^, and it has been 
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shown that the observed baryon asymmetry (cf. |[42||), 



B 



(0.6 - 1) ■ 10" 



(24) 



can be obtained for a wide range of parameters. 

In unified theories with right-handed neutrinos S — L is in general spontaneously 
broken. Unification also restricts the new parameters which are introduced by adding 
right-handed neutrinos to the standard model. In SO(IO) unification it is natural 
to assume a similar pattern of mixings and masses for leptons and quarks. This 
ansatz, together with the requirement of baryogenesis, also restricts the scale oi B — L 
breaking. The following discussion is closely related to |^|. 

The most general lagrangian for couplings and masses of charged leptons and 
neutrinos is given by 

- (25) 



Y 



-lL(pgieR - lL(pguJ^B. - ^^rM Ur + h.c. , 



where II = {i^l, ^l) is the left-handed lepton doublet and = (^9°, ip~) is the standard 
model Higgs doublet. The vacuum expectation value of the Higgs field (0) = w 7^ 
generates Dirac masses mi and m^, for charged leptons and neutrinos, 



mi = giv and mr, = QuV 



(26) 



which are assumed to be much smaller than the Majorana masses M. This yields 
light and heavy neutrinos 



with masses 



m^, ^ -K^'mo^m'^K* 



N 



mjsi ~ M . 



(27) 



(2^ 



as mass eigenstates. Here i^' is a unitary matrix which relates weak and mass eigen- 
states. Since the heavy neutrinos Ni are Majorana fermions, their decay to lepton 
and Higgs scalar violates lepton number. In the rest system the decay width of Ni 
reads at tree level, 



Di 



M ■ ^^^^ 

Interference between tree level and one-loop amplitudes (cf. Fig. 8) yields the CP 



{n' 0t + /) + r,, (at^ + 1) 



asymmetry |4l 



with /(x) = y/x 



1 - (1 + x) In 



J I] 

l + x 



M 



(30) 



X 
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The corresponding maximal B — L asymmetry is ei/ g*, where g* is the number of 



relativistic degrees of freedom (cf. W^ ) 



In order to generate the observed baryon asymmetry several conditions have to 
be fulfilled. First, the CP asymmetry ei has to be large enough; second, the out-of- 
equilibrium condition F/jj < nH{T = Mi) for the decaying heavy neutrino has to be 
fulfilled, where H is the Hubble parameter; third, the L violating interactions have 
to be sufficiently weak in order not to erase the generated lepton asymmetry. These 
conditions tend to favour small masses for the light neutrinos and a large scale of 
B — L breaking However, these constraints are not model independent. The 

corresponding bounds on the light neutrino masses can be considerably relaxed if 
appropriate chiral symmetries are effectively conserved in some temperature range in 
the symmetric phase [Q. 

All these conditions are automatically taken into account if one integrates the 
Boltzmann equations including all relevant interactions for the model under consider- 
ation. The results discussed below are based on such an analysis using the Boltzmann 
equations described in PT|. All three heavy neutrino families are taken into account 



as intermediate states whereas only the asymmetry generated by the lightest of the 
right-handed neutrinos is relevant, since the asymmetries generated by the heavier 
neutrinos are washed out. 



Neutrino masses and mixings 



Let us now consider a similar pattern of mixings and mass ratios for leptons and 
quarks, which is natural in SO(IO) unification. Such an ansatz is most transparent 
in a basis where all mass matrices are maximally diagonal. In addition to real mass 
eigenvalues two mixing matrices appear. One can always choose a basis for the lepton 
fields such that the mass matrices mi for the charged leptons and M for the heavy 




Figure 8: Contributions to the decay of a heavy Majorana neutrino 
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Majorana neutrinos A^, are diagonal with real and positive eigenvalues, 



mi 



/ me 


V 





m 








mr 



M 




(31) 



In this basis mo is a general complex matrix, which can be diagonalized by a biunitary 
transformation. Therefore, we can write m/5 in the form 



mn 




(32) 



where V and U are unitary matrices and the mj are real and positive. In the absence 
of a Majorana mass term V and U would correspond to Kobayashi-Maskawa type 
mixing matrices of left- and right-handed charged currents, respectively. 

According to Eqs. (|29| ) and ( pOD the CP asymmetry is determined by the mixings 
and phases present in the product m^m^), where the matrix V drops out. Therefore, 
to leading order, the mixings and phases which are responsible for baryogenesis are 
entirely determined by the matrix U. Correspondingly, the mixing matrix K in the 
leptonic charged current, which determines CP violation and mixings of the light 
leptons, depends on mass ratios and mixing angles and phases of U and V. Hence, 
there is no direct connection between the CP violation and generation mixing at high 
and low energies. 

Consider now the mixing matrix U. One can factor out five phases, which yields 



(33) 



where the Aj are the Gell-Mann matrices. The remaining matrix Ui depends on three 
mixing angles and one phase, like the Kobayashi-Maskawa matrix for quarks. In 



analogy to the quark mixing matrix we choose the Wolfenstein parametrization [§| 
as ansatz for Ui, 



I 



1-^ 

2 

-A 



A 



2 



\ A\^{\- p-ir]) -A\^ 



1 ) 



(34) 



where A and \p + ir]\ are of order one, while the mixing parameter A is assumed to 
be small. For the masses mi and Mj we assume a hierarchy like for up-type quarks. 



mi 
Ml = 



= hX^m^ 



m2 
M2 ■■ 



■- cX^m^j 
CX^Ms 



6,c = 0(1) 
B,C = 0{1) . 



(35) 
(36) 



For the eigenvalues m^ of the Dirac mass matrix this choice is motivated by SO (10) 
unification. The masses Mj cannot be degenerate, because in this case there exists a 
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basis for such that [7 = 1, which imphes that no baryon asymmetry is generated. 
For simphcity the masses Mj are assumed to scale hke the Dirac neutrino masses. 

The hght neutrino masses are given by the seesaw formula (^8|). The matrix 
K, which diagonalises the neutrino mass matrix, can be evaluated in powers of A. 
A straightforward calculation gives the following masses for the light neutrino mass 
eigenstates 



m,^ + O (A^) 



\C + e^'" B 
c2|C + e^^° B 



Ms ^ 



(37) 

(38) 
(39) 



The CP-asymmetry in the decay of the lightest right-handed neutrino Ni is easily 
obtained from Eqs. (BUp and (|5|)-(feB|), 



1 



BA^ 



A^^Im 



IGvr + \p + irj\'^ v 
This yields for the magnitude of the CP asymmetry, 

1 



(p _ ,^)2gi2(a+V3/3) 



C A 



(40) 



|ei| < 



BA^ 1 


P + irj\ 


2 


C2 + 


\P + i^l 


|2 



A^^ + 0(X' 



(41) 



How close the value of |ei| is to this upper bound depends on the phases a, j3 and 
Big{p + irj). Since ei oc m|/f^, one can already conclude that a large value of the 
Yukawa coupling m^/v will be preferred by this mechanism of baryogenesis. This 
holds irrespective of the neutrino mixings. 



Numerical results 



To obtain a numerical value for the produced baryon asymmetry, one has to specify 
the free parameters in the ansatz (^)-([36|). In the following we will use as a constraint 
the value for the z/^-mass which is preferred by the MSW explanation [^] of the solar 
neutrino deficit (cf. ||48[|), 

m^^ ~ 3 ■ 10"^ eV . (42) 

A generic choice for the free parameters is to take all 0{1) parameters equal to one 
and to fix A to a value which is of the same order as the A parameter of the quark 
mixing matrix, 

A = B = C = h = c=\p + ir]\r:^l, A ~ 0.1 . (43) 
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z = M,/- z = M,/- 

Figure 9: Time evolution of the neutrino number density and the B — L asymmetry 
for A = 0.1 and for = rrit (a) or = (h). The equilibrium distribution for 
Ni is represented by a dashed line, while the hatched area shows the measured value 
for the asymmetry. 



From Eqs. (§7p-(§g), (g^) and one now obtains, 

m^^ ~ 8 ■ 10"^ eV , m^^ ^ 0.15 eV . (44) 

Finally, a second mass scale has to be specified. In unified theories based on SO (10) 
the Dirac neutrino mass ms is naturally equal to the top-quark mass, 

mg = ^ 174 GeV . (45) 

This determines the masses of the heavy Majorana neutrinos Ni, 

M3 ~ 2 ■ 10^^ GeV , (46) 

and, consequently. Mi ~ 2 ■ 10^° GeV and M2 ~ 2 ■ 10^^ GeV. From Eq. (|T|) one 
obtains the CP asymmetry |ei| ~ 10~^, where we have assumed maximal phases. 
The solution of the set of Boltzmann equations discribed in PT| now yields the B — L 
asymmetry (see Fig. 9a), 

Yb-l ^ 3 ■ 10-1° , (47) 

which is indeed the correct order of magnitude. The precise value depends on un- 
known phases. 

The large mass M3 of the heavy Majorana neutrino A^3 (cf. (^6])), suggests that 
B — L is already broken at the unification scale Aqut ~ 10^^ GeV, without any 
intermediate scale of symmetry breaking. This large value of M3 is a consequence of 
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the choice ms ~ mj. To test the sensitivity of the result for Yb-l on this assumption, 
consider the alternative choice, rris = rrify ~ 4.5 GeV, with all other parameters 
remaining unchanged. In this case one obtains M3 = 10^^ GeV and |ei| = 5 ■ 10~^° 
for the mass of N3 and the CP asymmetry, respectively. Since the maximal B — L 
asymmetry is —ei/g* (cf. it is clear that the generated asymmetry will be too 

small. The solutions of the Boltzmann equations are shown in Fig. 9b. The generated 
asymmetry, Yb-l — 2 • 10~^^, is too small by more than two orders of magnitude. We 
conclude that high values for both masses and M3 are preferred, which is natural 
in SO(IO) unification. 

Models for dark matter involving massive neutrinos favour a r-neutrino mass 



ra^^ ~ 5 eV which is significantly larger than the value given in (^41). Such a 
large value for the r-neutrino mass can be accomodated within the ansatz described 
in this section. However, it does not correspond to the simplest choice of parameters 
and requires some fine-tuning. For the mass of the heaviest Majorana neutrino one 
obtains in this case M3 ~ 6 ■ 10^^ GeV. 

Without an intermediate scale of symmetry breaking, the unification of gauge 
couplings appears to require low-energy supersymmetry. This provides further sources 
for generating a B — L asymmetry [^, whose size depends on additional assumptions. 
In this case, especially constraints on the reheating temperature |^2[ and the possible 
role of preheating require further studies. 



4 Conclusions 

The observation that baryon and lepton number violating processes are in thermal 
equilibrium in the high-temperature phase of the standard model, is of crucial im- 
portance for the theory of baryogenesis. In particular it implies that the presently 
observed cosmological baryon asymmetry has been finally determined at the elec- 
troweak phase transition. 

During the past three years a quantitative understanding of the first-order elec- 
troweak phase transition for Higgs boson masses up to ttih ~ 70 GeV has been 
achieved by means of analytical and numerical methods. The transition to a crossover 
near mw ^ and the full understanding of the high-temperature phase still require 
further work. However, already now we know that for Higgs boson masses above the 
lower bound obtained at LEP, baryogenesis at the weak electroweak transition is very 
unlikely. 

Searching for alternatives to electroweak baryogenesis the connection between 
baryon number and lepton number in the symmetric phase is again crucial. It 
allows to generate the baryon asymmetry from a lepton asymmetry, as suggested 
by Fukugita and Yanagida. Necessary ingredients are right-handed neutrinos and 
Majorana masses, which appear naturally in SO(IO) unified theories. The exam- 
ple described in the previous section demonstrates that this mechanism can explain 
the observed cosmological baryon asymmetry without any fine-tuning of parameters. 
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Further sources for generating a lepton asymmetry exist in models with low-energy 
supersymmetry. All this suggests that, hoping for further progress in theory and new 
experimental results on neutrino properties, we can look forward to an intriguing 
interplay between non-pcrturbative processes in the standard model, early universe 
cosmology and neutrino physics. 
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